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The  physics  of  interacting  integer- spin  chains  has  been  a  topic  of  intense  theoretical  interest,  particularly 
in  the  context  of  symmetry-protected  topological  phases.  However,  there  has  not  been  a  controllable  model 
system  to  study  this  physics  experimentally.  We  demonstrate  how  spin-dependent  forces  on  trapped  ions 
can  be  used  to  engineer  an  effective  system  of  interacting  spin-1  particles.  Our  system  evolves  coherently 
under  an  applied  spin-1  XY  Hamiltonian  with  tunable,  long-range  couplings,  and  all  three  quantum  levels  at 
each  site  participate  in  the  dynamics.  We  observe  the  time  evolution  of  the  system  and  verify  its  coherence 
by  entangling  a  pair  of  effective  three-level  particles  (“qutrits”)  with  86%  fidelity.  By  adiabatically  ramping 
a  global  field,  we  produce  ground  states  of  the  XY  model,  and  we  demonstrate  an  instance  where  the  ground 
state  cannot  be  created  without  breaking  the  same  symmetries  that  protect  the  topological  Haldane  phase. 
This  experimental  platform  enables  future  studies  of  symmetry -protected  order  in  spin-1  systems  and  their 
use  in  quantum  applications. 
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I.  INTRODUCTION 

A  major  area  of  current  research  is  devoted  to  developing 
experimentally  controllable  systems  that  can  be  used  for 
quantum  computation,  quantum  communication,  and  quan¬ 
tum  simulation  of  many-body  physics.  To  date,  most 
experiments  have  focused  on  the  use  of  two-level  systems 
(“qubits”)  for  computation  and  communication  [1,2]  and 
for  the  study  of  spin- 1/2  (or  spinless)  many-body  phenom¬ 
ena  [3,4].  However,  there  are  a  variety  of  motivations  for 
performing  experiments  in  higher-dimensional  Hilbert 
spaces.  Contrary  to  the  intuition  that  enlarging  the  spin 
degree  simplifies  calculations  by  making  them  semiclass- 
ical  [5],  spin  >1/2  systems  inherently  have  more  com¬ 
plexity  and  cost  exponentially  more  resources  to  classically 
simulate.  For  instance,  it  is  computationally  easy  to  find  the 
ground-state  energy  of  a  spin- 1/2  chain  with  nearest- 
neighbor-only  interactions  in  one  dimension;  for  systems 
with  spin  7/2  or  higher,  the  problem  is  known  to  belong  to 
the  QMA-complete  complexity  class,  which  is  a  quantum 
analog  of  the  classical  NP-complete  class  [6,7].  The 
difficulty  of  this  problem  for  intermediate  spin  values, 
such  as  spin  1,  is  still  an  open  question.  From  a  more 
practical  point  of  view,  controllable  three-level  systems 
(“qutrits”)  are  useful  for  quantum  logic,  since  they  can 
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substantially  simplify  certain  operations  within  quantum 
algorithms  [8,9]  and  can  enhance  the  efficiency  of  quantum 
communication  protocols  [10]. 

When  individual  three-level  systems  are  coupled 
together,  they  can  be  used  to  encode  the  physics  of 
interacting  spin-1  particles.  Such  systems  have  attracted 
a  great  deal  of  theoretical  interest  following  Haldane’s 
conjecture  that  antiferromagnetic  Heisenberg  spin-1  chains, 
as  opposed  to  spin- 1/2  systems,  have  a  finite  energy  gap 
that  corresponds  to  exponentially  decaying  correlation 
functions  [11,12].  This  so-called  Haldane  phase  possesses 
a  doubly  degenerate  entanglement  spectrum  [13]  and  a 
nonlocal  string  order  [14,15],  which  is  related  to  the  order 
appearing  in  spin  liquids  [16]  and  in  the  fractional  quantum 
Hall  effect.  These  characteristics  suggest  that  the  Haldane 
phase  is  one  of  the  simplest  known  examples  of  a 
symmetry-protected  topological  phase  of  matter  [14]. 

In  addition  to  their  interesting  many-body  properties, 
topological  phases  may  be  exploited  in  a  more  applied 
setting.  The  Haldane  phase  is  useful  for  quantum  oper¬ 
ations  (for  instance,  as  a  perfect  quantum  wire)  [17,18]  and 
can  be  destroyed  only  by  crossing  a  phase  transition.  The 
finite  energy  gap  in  topological  spin-1  systems  makes  them 
a  potential  candidate  for  long-lived,  robust  quantum 
memories  [19],  and  schemes  using  symmetry-protected 
spin-1  phases  for  measurement-based  quantum  computa¬ 
tion  have  also  been  proposed  [20,21]. 

Several  groups  have  developed  controllable  three-level 
quantum  systems  by  using  pairs  of  photons  [22],  super¬ 
conducting  circuits  [23],  or  dressed  states  in  trapped  ions 
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[24]  to  implement  qutrits,  or  by  using  spinor  Bose-Einstein 
condensates  (BECs)  to  study  quantum  magnetism  [25-28]. 
However,  no  platform  has  yet  used  multiple  interacting 
qutrits  for  quantum  information  protocols  or  for  simulating 
lattice  spin  models.  In  this  paper,  we  use  trapped  atomic  ions 
to  simulate  a  chain  of  spin-1  particles  with  tunable,  long- 
range  XY  interactions  [29].  Our  system  performs  the  same 
basic  tasks  that  are  commonly  used  in  spin- 1/2  quantum 
simulations  [4,3CM-6],  such  as  observing  dynamical  state 
evolution  [44,45],  measuring  coherence  and  certifying  entan¬ 
glement  [31,40],  and  adiabatically  preparing  nontrivial 
ground  states  [42].  With  two  spin-1  particles,  we  observe 
coherent  evolution  under  the  XY  interactions  among  states  in 
a  “decoherence-free”  subspace  [47,48].  For  certain  states 
generated  by  the  XY  Hamiltonian,  we  can  verify  entangle¬ 
ment  between  a  pair  of  three-level  systems  with  fidelities  of 
up  to  86%.  Adding  a  time-dependent  global  field  allows  us  to 
adiabatically  prepare  the  ground  state  of  the  XY  model  for 
even  numbers  of  spins,  as  we  demonstrate  with  two-spin  and 
four-spin  chains.  For  odd  numbers  of  spins,  producing  the 
calculated  ground  state  is  not  possible  with  a  simple  adiabatic 
ramp  since  it  requires  crossing  a  first-order  phase  transition, 
hinting  at  the  existence  of  a  symmetry-protected  phase.  We 
observe  with  a  three-spin  chain  that  this  protection  of  the 
ground  state  is  preserved  even  in  the  presence  of  experimental 
imperfections.  It  is  straightforward  to  extend  our  techniques 
to  larger  collections  of  integer  spins,  providing  a  route  to 
accessing  higher-dimensional  Hilbert  spaces  experimentally. 
Additionally,  the  tools  demonstrated  here  could  enable  future 
studies  of  symmetry-protected  order  and  can  be  extended  to 
SU(3)  models  and  other  systems  of  higher  symmetry  [49]. 

II.  EXPERIMENTAL  IMPLEMENTATION 

The  spin-1  chain  is  represented  by  a  string  of  171Yb+ 
atoms  held  in  a  linear  Paul  trap.  Three  hyperfine  levels  in 
the  2 Si/ 2  ground  manifold  of  each  atom  are  used  to  encode 
the  spin-1  states:  |+)  =  \F  =  l,mF  =  1),  |-)  =  \F  =  1, 
mF  =  -1),  and  |0)  =  \F  =  0,  mF  =  0),  with  frequency 
splittings  of  oo±  between  the  |0)  and  |±)  states,  as  shown 
in  Fig.  1.  Here,  |+),  |— ),  and  |0)  are  the  eigenstates  of  Sz 
with  eigenvalues  +1,-1,  and  0,  respectively;  F  and  mF  are 
quantum  numbers  associated  with  the  total  angular  momen¬ 
tum  of  the  atom  and  its  projection  along  the  quantization 
axis,  defined  by  a  magnetic  field  of  ~5  G. 

We  apply  global  laser  beams  to  the  ion  chain  with  a  wave 
vector  difference  along  a  principal  axis  of  transverse 
motion,  polarized  to  drive  stimulated  Raman  transitions 
between  the  |0)  and  |-)  states  and  between  the  |0)  and  |+) 
states  with  equal  Rabi  frequencies  on  ion  i,  as  shown  in 
Fig.  1(a)  and  Fig.  1(b)  [50].  To  generate  spin-1  XY 
interactions,  we  apply  two  beat  frequencies  at  co_  +  fi 
and  coF  -  fi  to  these  respective  transitions,  where 
fi  is  slightly  detuned  from  the  transverse  motional 
frequencies,  as  shown  in  Fig.  1(c).  Under  the  rotating 
wave  approximations  oo±  » /u  and  within  the 
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FIG.  1.  (a)  Level  diagram  for  171  Yb+,  highlighting  relevant 

states,  (b)  Sketch  of  experimental  geometry,  showing  the  direc¬ 
tions  of  the  laser  wave  vectors  and  the  real  magnetic  field  relative 
to  the  ion  chain.  Both  beams  are  linearly  polarized,  one  along  the 
B  field  (providing  n  light)  and  one  orthogonal  to  the  B  field 
(providing  an  equal  superposition  of  o+  and  o~  light).  Multiple 
beat  notes  are  applied  by  imprinting  multiple  frequencies  onto 
one  beam  (in  this  case,  the  ^--polarized  beam),  (c)  Detailed  level 
diagram  of  the  251/2  ground  state,  showing  Raman  beat  notes  in 
relation  to  Zeeman  splittings  and  motional  sidebands  for  the 
center-of-mass  mode.  Level  splittings  are  not  to  scale. 

Lamb-Dicke  regime  (A k{xt)  1,  with  A k  the  wave  vector 
difference  of  the  Raman  beams  and  xt  the  position  operator 
of  the  ith  ion),  the  resulting  interaction  Hamiltonian  (with 
h  =  1)  is 

N  •  Q 

H  =  ini'mr  -  (-Si+amei(J,-m’n'>‘  +  SLale-^-^’). 
i,m=  1  ^v2 

(i) 

Here,  am  and  ah  are  the  phonon  operators  of  the  normal 
mode  m  with  frequency  com,  rji  m  =  bi  m  ^ h(Ak)2/ (2 Mcom) 
is  the  Lamb-Dicke  factor  (where  bi  m  is  the  normal  mode 
transformation  matrix  [51]  and  M  is  the  mass  of  a  single 
ion),  and  the  spin  raising  and  lowering  operators  Sl±  satisfy 
the  commutation  relations  [S+,Si]  =  2 Slz5ij.  In  the  limit 
where  the  beat  notes  are  far  detuned  ( \/a  -  com\) 
and  the  phonons  are  only  virtually  excited,  this  results  in  an 
effective  Hamiltonian  with  XF-type  spin-spin  interactions 
and  spin-phonon  couplings: 

Heif  =  J2J-f(s+sj-  +  sLsi) 

i<j 

+  J2ViJ(^m  +  l)Siz-(Siz)2].  (2) 

i,m 

The  pure  spin-spin  interaction  in  the  first  term  of  Eq.  (2) 
follows  the  same  formula  as  for  generating  spin- 1/2  Ising 
interactions  [50]: 
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hj  = 


m 


Vi,mH  j,m 

2(H  -  com) ' 


(3) 


When  i u  is  larger  than  the  transverse  center-of-mass 
frequency,  JLj  falls  off  with  distance  as  roughly 
J ij  ~  J()/\i  -  j\a,  where  J0  is  of  order  #1  kHz  and  a 
can  be  tuned  between  0  and  3  using  trap  and  laser 
parameters  [41,52]. 

The  V i  m  term  in  Eq.  (2)  is  given  by  a  similar  formula: 


(?1  i,m^i) 

8 (/<  -  a>m) ' 


(4) 


For  very  long-ranged  spin-spin  interactions  ( a  <  0.5),  or 
for  small  numbers  of  ions,  the  V i  m  terms  are  approximately 
uniform  across  the  spin  chain.  In  these  instances,  the  V i  m 
coefficient  can  be  factored  out  of  the  sum  over  ions  in 
Eq.  (2),  leaving  only  global  S\  and  (^)2  effective  field 
terms,  where  the  magnitude  of  the  S[  field  is  dependent  on 
the  number  of  phonons  in  each  mode.  The  phonon 
dependence  can  be  ignored  by  restricting  the  dynamics 
to  the  J2iS[  —  0  subspace,  as  discussed  below.  For  shorter- 
range  interactions  or  for  longer  chain  lengths,  the  V i  m 
terms  can  be  eliminated  completely  by  adding  an  additional 
set  of  beat  frequencies  at  co_  —  fi  and  oo+  +  //,  which  would 
generate  Ising-type  interactions  between  effective  spin-1 
particles  using  the  Mplmer-Sprensen  gate  [53].  As  has 
been  demonstrated  in  spin- 1/2  systems,  such  Ising  inter¬ 
actions  (along  with  a  strong  effective  magnetic  field)  can 
themselves  be  utilized  to  study  the  dynamics  of  XY 
interactions  [44,45]. 

The  ions  are  initialized  before  each  experiment  by 
cooling  the  transverse  modes  near  their  ground  state  of 
motion  ( n  «  0.05)  and  optically  pumping  the  spins  to  the 
1 00...)  state.  After  applying  the  Hamiltonian  in  Eq.  (2)  for 
varying  lengths  of  time,  we  measure  the  population  of  the 
state  |0)  at  each  site  by  imaging  spin-dependent  fluores¬ 
cence  [54]  onto  an  intensified  CCD  camera  and  observing 
which  ions  are  “dark.”  Because  both  of  the  |d=)  states 
appear  “bright”  during  the  detection  process  and  are 
scattered  into  an  incoherent  mixture  of  the  |  F  =  1)  states, 
our  current  setup  does  not  allow  discrimination  among  all 
three  possible  spin  states  in  a  single  experiment.  However, 
we  can  measure  the  population  of  either  |+)  or  |-)  by 
repeating  the  experiment  and  applying  a  n  rotation  to  the 
appropriate  |0)  |d=)  transition  before  the  fluorescence 

imaging.  For  instance,  measuring  an  ion  in  the  dark  state 
after  a  n  pulse  between  |0)  |+)  indicates  that  the  spin 

was  in  the  |+)  state  before  detection.  This  binary  discrimi¬ 
nation  is  not  a  fundamental  limit  to  future  experiments, 
since  populations  could  be  “shelved”  into  atomic  states  that 
do  not  participate  in  the  detection  cycle. 

Since  the  ions  are  initialized  to  the  1 00 . . . )  state,  and 
because  the  spin-spin  interactions  in  Eq.  (2)  conserve  the 
quantity  Sz,  the  dynamics  are  restricted  to  the  set 


of  states  with  Sz  =  0.  The  Sz  =  0  subspace  is  protected 
against  fluctuations  in  the  real  magnetic  field  A B(t),  which 
would  otherwise  result  in  an  unwanted  noise  term 
piBAB(t)Sz  (where  fiB  is  the  Bohr  magneton).  For  instance, 
the  T2  coherence  times  of  the  |0)  |d=)  transitions  are 

measured  to  be  0.5  ms,  limited  by  magnetic  field  noise. 
Nevertheless,  the  data  in  Figs.  2  and  3  exhibit  coherence 
and  entanglement  for  several  milliseconds  (limited  by  laser 
intensity  noise),  demonstrating  the  robustness  of  this 
decoherence-protected  subspace  against  time-varying 
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FIG.  2.  Dynamics  of  two  spin-1  particles  evolving  under  the 
XY  Hamiltonian  in  Eq.  (2).  (a)  We  measure  the  populations  for 
each  ion  to  be  in  the  state  |0) .  The  probability  of  both  ions  to  be  in 
|0)  (black  diamonds),  only  the  left  ion  in  |0)  (blue  triangles),  only 
the  right  ion  in  |0)  (green  squares),  and  neither  ion  in  |0)  (red 
circles)  are  plotted  together.  Similar  plots  for  the  |— )  state  (b)  and 
|+)  state  (c)  are  also  shown,  with  a  different  symbol  order 
emphasizing  that  the  red  circles  represent  the  same  state  in 
all  three  graphs.  The  dynamics  resemble  Rabi  flopping 
between  the  state  |00)  and  the  symmetric  superposition 
(|  4 — )  +  |  — \-))/y/2.  Solid  lines  represent  theoretical  predic¬ 
tions,  with  the  only  free  parameters  being  the  magnitude  of  the  Sz 
and  (Sz)2  gradients  discussed  in  the  text.  In  (c),  the  interaction 
Jx  2  drifts  at  long  times  compared  to  that  estimated  from  Eq.  (3), 
consistent  with  an  observed  drift  in  radial  trap  frequencies  during 
the  data  collection.  Error  bars  (ls.d.)  show  the  statistical 
uncertainty  based  on  500  repetitions  of  the  experiment. 
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FIG.  3.  (a)  Illustration  of  the  rotations  performed  before 

measuring  the  parity,  where  rotations  by  <9  and  cp  are  defined 
in  Eq.  (5).  Also  shown  are  the  ideal  initial  state  and  the  states 
produced  after  each  step,  (b)  The  parity  of  the  final  state  oscillates 
as  a  function  of  the  final  pulse  phase  cp.  Fitting  the  function 
C  —  A  cos  2 cp  +  B  sin  2 cp  (red  dashed  line)  to  the  data  results  in  an 
amplitude  A  =  0.86  >  0.5,  demonstrating  entanglement,  (c)  Am¬ 
plitude  of  the  parity  oscillation  after  various  durations  of  the  spin- 
spin  interaction,  showing  the  peak  amplitude  for  each  time 
~(2 n  +  1) / (2y/2J i  2) .  The  data  in  (b)  correspond  to  the  high- 
est-contrast  point  in  (c).  The  dashed  line  is  a  guide  to  the  eye, 
suggesting  the  expected  behavior  of  sinusoidal  oscillations 
between  a  product  state  and  an  entangled  state,  along  with  decay 
due  to  decoherence  at  longer  times.  The  chosen  durations  are  not 
evenly  spaced  past  0.6  ms:  for  each  point,  the  duration  is  chosen 
such  that  the  population  in  |00)  is  at  a  local  minimum,  and  drifts 
in  the  radial  trap  frequencies  lead  to  small  changes  in  Jl2. 

magnetic  fields.  We  expect  that,  independent  of  laser- 
induced  decoherence,  the  coherence  times  within  this 
subspace  could  reach  seconds  [54],  since  the  magnetic 
field  noise  A B(t)  will  introduce  the  same  second-order 
effects  as  in  the  case  of  the  hyperfine  “clock”  qubit.  This 
improved  coherence  does  not  require  any  extra  overhead 
for  either  initialization,  readout,  or  generation  of  spin-spin 
interactions.  Additionally,  remaining  within  the  Sz  =  0 
subspace  does  not  substantially  limit  the  size  of  the 
accessible  Hilbert  space,  since  the  number  of  states  in 
the  Sz  =  0  subspace  of  N  spin-1  particles  scales  as 


~3N /  {2y/N)  for  large  N,  which  is  exponentially  greater 
than  the  2N  states  accessible  in  a  spin- 1/2  system. 

Our  implementation  of  a  decoherence-protected  subspace 
is  closely  analogous  to  techniques  that  have  been  used  for 
extending  the  coherence  time  of  certain  spin- 1/2  quantum 
simulations  [45],  and  for  creating  a  logical  qubit  that  is 
robust  to  environmental  perturbations  [48],  for  which  a 
universal  gate  set  has  been  demonstrated  [55].  In  addition  to 
the  quantum  simulations  discussed  in  this  work,  our 
approach  could  potentially  be  used  in  quantum  logic  to 
encode  a  logical  qutrit  using  two  physical  qutrits.  For 
example,  we  would  encode  |+)L  =  ||j  4 — ),  |0)L  =  |00), 
and  |-)L  =  |  —  +).  Such  a  logical  qutrit  could  be  trivially 
initialized,  deterministically  measured  in  a  single  shot, 
without  recourse  to  shelving  techniques,  by  measuring 
the  |+)  or  |— )  components  of  the  physical  qutrits  (using 
n  rotations  from  |d=)  to  |0)),  and  operated  upon  using  the 
AT- type  interactions  demonstrated  here  (which  would  be 
equivalent  to  a  logical  Sx  or  Sy ,  depending  on  laser  phases) 
or  a  local  Sz  applied  to  one  of  the  physical  qutrits 
(equivalent  to  a  logical  Sz).  A  conditional  phase  gate 
between  two  logical  qutrits,  which  (along  with  single-qutrit 
manipulations)  permits  a  universal  gate  set  for  quantum 
computation  with  qutrits  [56],  could  be  accomplished  by  the 
application  of  an  Ising-type  (, SlzS{ )  interaction  involving  one 
physical  qutrit  from  each  pair.  These  operations  could  be 
accomplished  on  similar  time  scales  to  those  for  typical  two- 
qubit  gates  in  trapped  ions,  on  the  order  of  hundreds  of  jus. 
We  note  that  a  complementary  technique  has  been  demon¬ 
strated  for  encoding  a  decoherence-resistant  qutrit  (or  qubit) 
in  a  single  atom,  at  the  expense  of  requiring  the  continuous 
application  of  microwave  dressing  fields  [24,57]. 

III.  COHERENT  DYNAMICS  OF  TWO  SPINS  AND 
ENTANGLEMENT  VERIFICATION 

For  a  system  of  two  spins,  dynamical  evolution  under  the 
Hamiltonian  in  Eq.  (2)  can  be  understood  as  Rabi  flopping 
between  the  |00)  and  (|  4 — )  +  |  — h))/v^  states  with 
Rabi  frequency  y/2 Ji2.  This  behavior  is  shown  in  Fig.  2, 
where  Figs.  2(a)-2(c)  show  the  probability  of  each  ion  to  be 
in  the  |0),  |— ),  and  |+)  states,  respectively.  The  population 
remains  in  the  Sz  =  0  subspace,  as  expected:  Fig.  2(a) 
shows  the  absence  of  the  Sz  ^  0  states  (|0+),  |0 — ),  |  +  0), 
and  |  -  0)),  while  Figs.  2(b)  and  2(c),  respectively,  show 

the  absence  of  the  other  Sz  ^  0  states  | - )  and  |  +  +). 

The  drift  in  JX  2  evidenced  in  Fig.  2(c)  could  be  stabilized  in 
future  experiments  by  feeding  back  to  the  trap  rf  voltage  to 
better  stabilize  the  radial  trap  frequencies. 

The  different  ions  i  can  experience  position-dependent 
S[  and  (^)2  shifts,  whose  effects  are  visible  in  the  slight 
divergence  of  the  blue  and  green  curves  in  Figs.  2(b) 
and  2(c).  We  attribute  this  effect  to  a  micromotion  gradient, 
since  the  shifts  can  be  compensated  by  adjustments  of 
the  voltages  on  the  dc  trap  electrodes,  and  since  (based  on 
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calibration  experiments  with  three  ions  in  which  we 
measure  the  transition  energies  from  |000)  to  other  states 
in  the  Sz  =  0  subspace  in  the  presence  of  an  uncompen¬ 
sated  gradient)  the  energy  levels  are  consistent  with  a  linear 
gradient  in  position. 

The  calculation  overlaid  in  Fig.  2  includes  the  site- 

dependent  terms  (200  Hz)*#  +  (150  Hz ,)(s£2))2,  which 
were  left  as  free-fitting  parameters  when  numerically 
evolving  the  Schrodinger  equation  under  the 
Hamiltonian  in  Eq.  (2),  since  the  calibrations  mentioned 
above  have  a  precision  of  only  ~200  Hz.  The  plotted 
curves  assume  strictly  unitary  evolution  (i.e.,  no 
decoherence)  over  the  time  scale  of  the  experiments. 

At  a  time  t  =  0.5/(a/27i,2)>  which  is  roughly  0.27  ms  in 
Fig.  2,  the  system  is  left  approximately  in  the  entangled 
state  ((  +  —)  +  [ — |-)) /V 2.  To  verify  entanglement  in 
the  system,  one  could  use  spin-1  analogs  of  Bell-type 
inequalities  [58],  which  require  many  local  rotations  but 
are  sensitive  to  maximally  entangled  states  like 
(| 00)  +  |  +  — )  +  |  —  +))/v/3*  However,  for  the  class  of 
states  generated  by  the  XY  interactions,  a  much  simpler 
series  of  global  rotations  is  sufficient  to  verify  entangle¬ 
ment.  The  analysis  consists  of  performing  three  sequential 
rotations  on  the  |0)  to  |=b)  transitions, 

R0±(e,  4>)  =  e[(^/2)Et[^(i±xoi),+^^(iox±i)d])  (5) 

before  measuring  the  population  in  |0).  The  rotation 
sequence  is  given  by  R<d+(n/2,  (-tt/2,  0)R0_(7i,  0), 

with  the  rotations  applied  from  right  to  left.  The  first 
two  rotations  map  the  state  (|  H — )  +  |  — b))/\/2  to 
(|00)  +  |  +  +))/\/2,  while  the  phase  of  the  third  rotation 
is  varied  to  analyze  the  entanglement  of  this  resulting  state 
[59].  The  parity  n  =  J2j=o(~^)jPj  (with  Pj  the  probability 
of  j  atoms  in  |0))  oscillates  as  a  function  of  the  phase  cp  of 
the  third  pulse,  and  the  amplitude  of  its  oscillation  depends 
on  the  off-diagonal  density  matrix  elements: 

n (cp)  =  C  -b  ^cos  2cp(P__  +  P++  -  P+_  -  P_+ 

-2|p+- -+|  -  2 |p — ,++|)  +~^sm.2(p{2\p_+>++\ 

+  2|/>+_ -h-I  -  2\p — _+|  -  2|p_+j — |),  (6) 

where  Pt  is  the  population  in  state  |/)  (|/)  =  j - ),  |  — b), 

etc.),  pi  j  is  the  off-diagonal  density  matrix  element 
quantifying  the  coherence  between  |/)  and  | j),  and  C  is 
a  constant  offset  that  depends  on  the  various  density  matrix 
elements  but  not  on  the  phase  cp  of  the  final  rotation.  The 

populations  in  |  +  +)  and  | - )  are  negligible,  simplifying 

this  expression: 

n (cp)  «  C  -  A  cos  2cp,  (7) 


where  the  oscillation  amplitude 

A  =  ^  {p+-  +  P-+  +  2l/U-,-+l)  (8) 

is  akin  to  the  entanglement  fidelity  T  of  Greenberger- 
Horne-Zeilinger  (GHZ)  states  in  two-level  systems  [59]. 
Measuring  the  amplitude  A  of  the  parity  oscillation  n [cp) 
then  allows  us  to  verify  entanglement  for  certain  classes  of 
states.  According  to  an  analysis  analogous  to  that  in 
Ref.  [59],  the  following  inequality  holds  for  all  separable 
qutrit  states: 

2A  +  P0o  +  2|p+- ool  +  2|/?-+,ool  ^  1-  (9) 

Hence,  violation  of  this  inequality  demonstrates  entangle¬ 
ment  between  spin-1  particles  or  qutrits,  and  measuring  an 
amplitude  of  A  >  1/2  is  sufficient  to  violate  the  inequality. 

Figure  3(b)  shows  an  example  of  the  measured 
parity  curve  used  to  extract  the  amplitude  A  and  verify 
entanglement  between  the  qutrit  pair.  Such  measurements 
can  be  repeated  for  different  durations  of  exposure  to  the 
XY  Hamiltonian.  At  times  t  =  (2  n  +  1)  /  (2\/2J  l2) 
(n  =  0, 1,2, ...),  the  system  should  again  be  in  the  state 
(|  H — )  +  |  — b))/ y/2,  while  at  times  t  =  n/{y/2J\^),  it 
should  return  to  the  unentangled  product  state  |00).  The 
result  is  plotted  in  Fig.  3(c). 

Two  known  sources  of  dephasing  contribute  to  the 
observed  loss  of  coherence  in  the  experiment.  First,  laser 
intensity  fluctuations  and  pointing  instability  cause  noise  in 
the  spin-spin  coupling  term,  leading  to  apparent  dephasing 
when  many  repetitions  are  averaged  together.  These 
fluctuations  could  be  compensated  in  future  experiments 
by  variants  of  the  method  of  composite  pulses  [60,61].  The 
second  dephasing  source  results  from  inhomogeneities 
in  the  V  im  term  [Eq.  (2)]  across  the  chain,  which  will 
cause  different  spins  to  acquire  phases  at  different  rates. 
This  is  mainly  a  concern  in  chains  longer  than  two  ions 
(although  inhomogeneous  Rabi  frequencies  £2*  can  cause 
V i  m- induced  dephasing  even  for  two  spins),  and  could  be 
compensated  by  adding  an  extra  driving  term  to  cancel  the 
inhomogeneities  or  by  applying  a  series  of  echo  pulses 
[62,63].  Fluctuating  external  magnetic  fields  and  off- 
resonant  coupling  to  the  carrier  transition  would  ordinarily 
add  dephasing  noise  along  the  z  direction,  but  these 
have  been  suppressed  here  by  working  in  the  Sz  =  0 
subspace. 

IY.  GROUND  STATE  PRODUCTION 

We  can  also  add  an  effective  (5/)2  field  term, 

(^l)2’  to  the  Hamiltonian  by  shifting  the  beat 
frequencies  of  the  Raman  lasers  to  cd+-  p-  D  and 
co_  p  —  D.  This  (S?,)2  term  can  be  used  to  adiabatically 
prepare  the  ground  state  of  the  XY  Hamiltonian  in  Eq.  (2). 
As  before,  the  spins  are  prepared  in  |00 _ ),  which  is  the 
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FIG.  4.  Measurements  of  the  prepared  two-spin  states  (narrow 
blue  bars)  after  ramping  an  ( Sz )2  field,  compared  to  the  values 
expected  for  the  calculated  ground  state  (gray  bars).  As  in 
Fig.  2(a)-(c)  show  the  measured  populations  when  the  dark  state 
is  set  to  be  |0),  |— ),  or  |+),  respectively. 


approximate  ground  state  of  Eq.  (2)  in  the  presence  of  a 
large  (5  kHz)  (S^)2  field.  This  field  is  then  ramped  down 
slowly  according  to  D(t)  =  (5  kHz)^-^0-167  ms).  Figures  4 
and  5  show  the  populations  measured  at  the  end  of  the  (Slz)2 
ramp  for  two  and  four  spins,  which  match  reasonably  well 
with  the  calculated  ground  state. 

Measurements  of  populations  in  the  Sz  basis  necessarily 
discard  phase  information  about  components  of  the  final 
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FIG.  5.  Measurements  of  the  prepared  four-spin  states  (narrow 
blue  bars)  after  ramping  an  ( Sz )2  field,  compared  to  the  values 
expected  for  the  calculated  ground  state  (gray  bars).  Again, 
(a)-(c)  show  the  measured  populations  when  the  dark  state  is  set 
to  be  |0),  |— ),  or  |+),  respectively.  Here,  the  interaction  range  is 
given  by  Juj  ~  l/|i  -  ;'|0-3. 


FIG.  6.  Following  an  adiabatic  ramp,  the  parity  of  the  final 
state  is  measured  as  a  function  of  the  final  rotation  phase  cp  (see 
text  for  rotation  protocol).  Dashed  and  dot-dashed  lines  represent 
the  theoretically  expected  values  for  the  ground  state, 
|00)/\/2  -  (|  -  +)  +  |  +  -))/2,  and  highest  excited  state, 
100) /x/2  -F  (|  — h)  +  |  H — ))/2,  respectively.  The  phase  of  the 
oscillation  reveals  that  the  relative  phases  in  the  prepared  state  are 
consistent  with  the  expected  ground  state. 


state.  This  can  be  important  in  many  spin  models,  including 
the  XY  model,  where  such  measurements  alone  cannot 
discriminate  between  different  eigenstates.  For  example, 
the  ground  state  of  an  XY  model  with  two  spin-1  particles  is 

|00)/\/2  -  (|  — h)  +  |  H — ))/2,  while  the  highest  excited 
state  is  |00) / \fl  +  (|  -  +)  +  |  +  -))/2,  differing  by  only  a 
relative  phase.  We  check  that  we  are  creating  the  ground 
state  after  our  adiabatic  protocol  by  applying  a  pair  of 
rotations,  R^n/2,  0),  and  measuring  the 

parity  n,  as  was  done  in  the  entanglement  analysis.  This 
is  expected  to  result  in  Yl{cp)  =  (3/8)  d=  (1/2)  cos  cp ,  where 
+  and  -  correspond  to  the  ground  and  highest  excited 
states,  respectively.  As  shown  in  Fig.  6,  our  measurements 
are  consistent  with  having  prepared  the  two-spin 
ground  state. 
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FIG.  7.  Measurements  when  the  |0)  state  is  set  dark  (a),  the  |+) 
state  is  dark  (b),  and  the  | — )  state  is  dark  (c),  of  the  prepared  three- 
spin  state  after  adiabatically  ramping  a  global  (S/)2  field  (narrow 
blue  bars).  The  data  agree  closely  with  the  calculated  populations 
in  the  first  excited  state  (gray  bars),  while  showing  little  overlap 
with  the  expected  populations  in  the  ground  state  (wide,  hatched 
red  bars). 
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V.  TOWARD  HALDANE  PHYSICS 

A  more  long-term  goal  for  spin-1  quantum  simulations  is 
to  produce  and  study  ground  states  in  the  Haldane  phase 
[29].  It  is  known  that  an  XY  model  with  both  nearest- 
neighbor  and  next-nearest-neighbor  interactions  can  exhibit 
a  symmetry-protected  Haldane  phase  [64],  and  it  remains 
an  open  question  whether  a  generic  long-range  XY  model 
would  show  the  same  behavior.  Already  with  our  exper¬ 
imentally  implemented  Hamiltonian,  we  find  a  useful  test 
case  where  the  symmetry  of  the  ground  state  prevents  it 
from  being  created  via  the  simple  adiabatic  protocol 
described  above. 

The  ground  state  |i j/)  of  a  long-range  XY  model  can  be 
calculated  exactly  for  three  spins.  For  our  experimental 
coupling  strengths  JLj  —  1/| i  —  j |036, 

k)gs  =  >/0.16(|0  — h)  —  |0  H — )  +  |  —  +0)  -  |  +  -0)) 

+  V0l8(|+0-)-|-0+)).  (10) 

This  state  has  a  99.9%  overlap  with  a  three-spin  Affleck- 
Kennedy-Lieb-Tasaki  (AKLT)  state  [65],  which  is  the 
canonical  example  of  a  ground  state  in  the  Haldane  phase 
that  can  be  written  in  closed  form  for  any  number  of  spins. 
The  state  in  Eq.  (10)  is  antisymmetric  with  respect  to  the 
same  symmetries  that  govern  the  Haldane  phase,  such  as 
left-right  spatial  inversion  of  the  chain  or  a  global  rotation 
about  Sx  by  n  (which  sends  |+)  to  |— )  and  vice  versa). 
However,  since  the  starting  state  |000)  and  the  applied 
Hamiltonian  are  symmetric  with  respect  to  these  opera¬ 
tions,  we  should  be  unable  to  reach  the  antisymmetric 
ground  state  with  a  simple  adiabatic  ramp.  Indeed,  we  find 
numerically  that  a  first-order  phase  transition  separates  the 
symmetric  and  antisymmetric  ground  states,  which  cannot 
be  adiabatically  connected  without  breaking  inversion  and 
rotational  symmetry.  For  the  three-spin  experiment  in 
Fig.  7,  we,  hence,  prepare  a  state  close  to  the  first  excited 
state  rather  than  the  ground  state.  This  observation  suggests 
that  even  in  the  presence  of  various  experimental  imper¬ 
fections,  the  ground  state  of  our  three-spin  XY  model 
enjoys  the  same  symmetry  protection  as  the  Haldane  phase. 

In  this  paper,  we  demonstrate  the  basic  ingredients  that 
are  needed  for  the  implementation  of  quantum  simulations 
with  spins  greater  than  1/2,  using  a  platform  in  which  there 
are  no  fundamental  limitations  to  scaling  toward  larger  spin 
chains.  We  believe  that  this  work  opens  paths  for  studying 
the  exciting  physics  beyond  spin- 1/2  systems,  and  we  have 
already  taken  the  first  steps  towards  exploring  the  richness 
of  topological  phases.  In  particular,  for  a  long-range  spin-1 
XY  model,  we  demonstrate  coherent  Schrodinger  evolution 
and  the  capability  to  create  symmetric  ground  states.  We 
observe  that  for  odd  numbers  of  spins,  symmetry  consid¬ 
erations  prevent  us  from  creating  ground  states  that  bear  a 
close  resemblance  to  AKLT  states  and,  hence,  may  belong 
to  the  Haldane  phase.  Future  work  will  address  the 


questions  of  how  to  add  a  Heisenberg  term  and  sym¬ 
metry-breaking  perturbations  to  the  Hamiltonian  so  as  to 
prepare  antisymmetric  ground  states  [29],  which  will  allow 
us  to  create  and  probe  interesting  edge  states  in  the 
Haldane  phase. 
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